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(Dated: February 7, 2008) 

We report the results of a search for T, CP and CPT violation in B^-B^ mixing using an inclusive 
dilepton sample collected by the BABAR experiment at the PEP-II B Factory. Using a sample of 232 
million BB pairs, with a simultaneous likelihood fit of the same-sign and opposite-sign dileptons, 
we measure the T and CP violation parameter \q/p\ — 1 = (—0.8 ± 2.7(stat.) ± 1.9(syst.)) x 10""^, 
and the CPT and CP parameters Imz = (-13.9 ± 7.3(stat.) ± 3.2(syst.)) x 10"^ and AF x Rez = 
( — 7.1 ± 3.9(stat.) ± 2.0(syst.)) x 10~^ps~^. The statistical correlation between the measurements 
of Imz and AF x Rez is 76 %. 

PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 



Since the first observation of CP violation in 1964 
the neutral kaon system has provided many results prob- 
ing the CPT and T discrete symmetries in K'^-K'^ 
mixing. Similarly, the BABAR experiment can investigate 
T, CP, and CPT violation in B^-B^ mixing. 

The physical states (solutions of the complex effective 
Hamiltonian for the B^-B^ system) can be written as 



\Bh) = px/miB") -9x/T^|B"). 

where H and L stand for Heavy and Light. In the case 
of CPT invariance, the complex parameter z is equal to 
0. Similarly, T invariance leads to \q/p\ = 1. Finally, CP 
invariance requires both \q/p\ = 1 and z = 0. 

Inclusive dilepton events, where both B mesons decay 
semileptonically b Xtv {I — e or fi), represent 4% of 
all T(45) BB decays and provide a very large sample 
to study T, CPT and CP violation in mixing. In the 
direct b ^ I decay process, the flavor B^{B^) is tagged 
by the charge of the lepton {l~ ) . 

At the T(4S') resonance, neutral B mesons are 
produced in a coherent p-wave state. At the instant 
that the first B meson decays, the second B meson has 
the opposite flavor. Then, the second B meson will 
continue to evolve in time. Defining the time difference 
as ISX — — t~ where t~^{t~) is the decay time of the 
neutral B tagged by £^ {£~ ) , and neglecting second order 
terms in z, the decay rates for the three configurations 
{£+£+, £-£- and £+£-) are given by 



++ 



N 



N~ 



N+- 



2 'g 



-r|At| 



f , , Al Zlt, , , ^ , -1 

|cosh( — - — ) — cos(AmAt) J-, 

,2 r , , AFAt, , ^ ^ , 1 

I I cosh ( — - — ) - cos(AmAt)|, 

AFAt, . , ,AFAt, 



I cosh ( 



2 L '2 

-I- cos 



2 Rez sinh(- 



(ArriAt) -I- 2Imz sin(AmAt)|, (1) 



where Am is the B'^-B'^ oscillation frequency, F is the 
average neutral B decay rate and AF is the decay rate 
difference between the two physical states. 

The same-sign dilepton asymmetry Arp/Qp, between 
the two oscillation probabilities P{B'^ B°) and 



P[B^ B^) probes both T and CP symmetries and 
can be expressed in terms of \q/p\'. 



A 



T/CP 



P{B^ B^) - P{B^ B°) 
P{BO BO) 

N++ -N— 



N++ + N- 



P{BO B°) 

1 - \q/p\^ 
" l + \qlp\^' 



(2) 



Standard Model calculations |j| predict the size of this 
asymmetry to be at or below 10^"^. A large measured 
value would be an indication of new physics. 

Similarly, the opposite-sign dilepton asymmetry, 
A(j prp , between events with At > and At < com- 
pares the B^ B^ and B'^ B" probabilities and is 
sensitive to CPT and CP violation. This asymmetry is 
given by 

P(B" ^ B") - P(B" -> W) 



A. 



CPT (CP 



(lAt|) 



P(BO PO) + P(P0 ^ PO) 
jV+~ (At > 0) - Af+~ (At < 0) 
iV+-(At > 0) -l-iV+-(At < 0) 
Imzsin(AmAt)-Rezsinh(^^) 



cosh(- 



+ cos(AmAt) 



.(3) 



As |AF|/F < 1 [g, we have Rezsinh(ArAt/2) ~ AF x 
Rez x (At/2) and this asymmetry is not sensitive to the 
CPT-violating term Rez alone, but to the product AF x 
Rez. 

In this Letter, we present measurements of \qlp\; Imz 
and AF x Rez with a simultaneous likelihood fit of the 
same-sign and opposite-sign dilepton At distributions. In 
the cosh(AFAt/2)term, we fix |AF| to 0.005 ps"\ the 
value reported in [2| with a 90% confidence-level limit of 
0.055 ps-i. 

This study is performed with events collected by the 
PABAR detector |5j at the PEP-II asymmetric-energy B 
Factory between October 1999 and July 2004. The in- 
tegrated luminosity of this sample is about 211 fb~^ 
recorded at the T(4S') resonance ( "on-resonance" ) (232 
million BB pairs) and about 16 fb~^ recorded about 40 
MeV below the T(4S') resonance ( "off- resonance" ) . 

The event selection is identical to that described in . 
Non-PP events are suppressed by applying requirements 
on the ratio of second to zeroth order Fox- Wolfram mo- 
ments 0, the squared invariant mass, the aplanarity and 
the number of charged tracks of the event. 
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Lepton candidate tracks must have at least 12 hits in 
the drift chamber, at least one ^-coordinate hit in the sili- 
con vertex tracker (SVT), and a momentum in the T{AS) 
center-of-mass system between 0.8 and 2.3 GeV/c. Elec- 
trons are selected by requirements on the ratio of the 
energy deposited in the electromagnetic calorimeter and 
the momentum measured in the drift chamber. Muons 
are identified through the energy released in the calorime- 
ter, as well as the strip multiplicity, track continuity, and 
penetration depth in the instrumented flux return. Lep- 
ton candidates are rejected if their signal in the detector 
of internally reflected Cherenkov light is consistent with 
that of a kaon or a proton. The electron and muon se- 
lection efficiencies are about 85% and 55%, with pion 
misidentification probabilities around 0.2% and 3%, re- 
spectively. 

Electrons from photon conversions are identified and 
rejected with a negligible loss of efficiency for signal 
events. Leptons from J/ijj and iJj{2S) decays are identi- 
fied by pairing them with other oppositely-charged can- 
didates of the same lepton species, selected with looser 
criteria. 

The separation between direct leptons (6 — > £) and 
background from the b ^ c ^ £ decay chain {cascade 
decays) is achieved with a neural network that combines 
five discriminating variables: the momenta and opening 
angles of the two lepton candidates, and the total visible 
energy and missing momentum of the event, all computed 
in the T{AS) rest frame. 

Of the original sample of 232 million BB pairs, 1.4 
million pass this dilepton selection. 

Since the asymmetry At/cp is expected to be small, 
we have determined the possible charge asymmetries in- 
duced by the detection and reconstruction of electrons 
and muons. The charge asymmetries are defined by 
a = e^)/(e+ + e^) where e~^{e~) is the efficiency for 
positive and negative particles. As the lepton efhciencies 
and purities depend on their allowed phase space, we con- 
sider separately the asymmetry for the higher and lower 
momentum lepton, respectively, oi and a2. 

The charge asymmetry of track reconstruction is mea- 
sured in the data by comparing tracks reconstructed us- 
ing only the SVT with those passing the dilepton track 
selection, obtaining atrk = (0.8 ± 0.2) x 10"'^. 

The lepton identification efficiencies are measured as a 
function of total momentum and polar and azimuthal an- 
gles, with a control sample of radiative Bhabha events for 
electrons, and with a,ee —^ fifij control sample for muons. 
The misidentification probabilities are determined with 
control samples of kaons produced in 7r+_D° 
7r+_R'^7r+ (and charge conjugate) decays, pions produced 
in Ks TT+TT" decays, one-prong and three-prong r de- 
cays, and protons produced in A decays. 

The control samples show that the muon track re- 
construction efficiency has a charge asymmetry reaching 
~ 5 X 10~^ and that the positive kaons are more likely 



than negative kaons to be misidentified as muons at the 
20-30% level. As a consequence, in the likelihood fit (de- 
scribed below), we float the charge asymmetries ajj""" and 
ijcasc (jij.g(2^ a^jjjj cascade muons. 

For electrons, the charge asymmetry averaged over the 
signal phase space is Ue = (0.4±0.2) x 10~^ and we find 
that antiprotons with momentum ~ 1 GeV/c are signifi- 
cantly more likely than protons to be misidentified, due 
to annihilation with nucleons in the calorimeter material. 
Based on the charge asymmetry in tracking and in identi- 
fication, we fix the charge asymmetry for the direct elec- 
trons with the higher momentum to a'^Y = 1-2 x 10""^. 
For the lower momentum direct electrons and the cas- 
cade electrons, for which antiprotons contamination is 
more important, we correct the initial charge asymmetry 
by the fraction of antiprotons estimated with generic BB 
Monte Carlo samples and the proton control sample, this 
gives the following charge asymmetries: af**" = 0.8x10^'^, 
a™"^ = 0.5 X 10-3, and a^J^'" = 0.2 x lO'^. 

In the inclusive approach used here, the z coordinate 
of the B decay point is approximated by the z position of 
the point of closest approach between the lepton candi- 
date and an estimate of the T(4S') decay point in the 
transverse plane. The T(45') decay point is obtained 
by fitting the two lepton tracks to a common vertex in 
the transverse plane that is constrained to be consistent 
with the beam-spot position. The proper time difference 
At between the two B meson decays is determined from 
Az = z^ — z~ , the difference in z between the leptons 
and by At = Az/{(3"f)c with a nominal Lorentz 
boost {f3j) — 0.55. In case of same-sign dileptons, the 
sign of At is chosen randomly. 

We model the contributions to our sample from BB de- 
cays using five categories of events, i, each represented by 
a probability density function (PDF) in At, V"''^. Their 
shapes are determined using the i3°B° (n) and B^B^ (c) 
Monte Carlo simulation separately, with the approach de- 
scribed in j^]. 

The five categories are the following. First, the pure 
signal events with two direct leptons (sig) represent 81% 
of the BB events and give information on the T, CPT 
and CP parameters. Then, we consider two categories of 
cascade decays: those with a direct lepton and cascade 
lepton from the opposite B decays (o6c), and those with 
direct lepton and cascade lepton from the same B decay 
(sbc). According to generic BB Monte Carlo simulation, 
their contributions are around 9% and 4% respectively. 
In addition, 3% of the dilepton events originate from the 
decay chain b ^ t~ £~ (Idlr) which tags the B flavor 
correctly. Finally, the remaining events {other) consist 
mainly of one direct lepton and one lepton from charmo- 
nium resonances in the B decays. 

The sig event PDF, V^-g, is obtained by the convolu- 
tion of an oscillatory term containing the T, CPT and 
CP parameters (Eq. ^ for neutral B decays or an expo- 
nential function for charged B decays, with a resolution 
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function which is the sum of three Gaussians. The widths 
of the core and tail Gaussians and the fractions of core 
and outher Gaussians are free parameters in the fit. The 
width of the outher Gaussian is fixed to 8 ps. The means 
of the Gaussians are fixed to zero |^ . 

The obc event PDF, V^^^, is modeled by the convo- p 
lution of At-dependent terms of a form similar to those 
of the signal with a resolution function which takes into 
account the effect of the charmed meson lifetimes. Since 
both short-lived and Ds, and long-lived D'^ mesons 
are involved in cascade decays, the resolution function 
for the long-lived and short-lived components is a sum of 
three Gaussians, which are convoluted with double-sided 
exponentials. To correct the effect of possible outliers not 
observed with the Monte Carlo simulation, the fraction 
of the third Gaussian is free in the fit. Similarly, we take 
the effect of the charmed mesons into account in the sbc 



event PDF, V] 



sbc ■ 



The PDF for Idlr events, V^f^^^ is similar to that of 
the sig events. The resolution function used takes into 
account the r lifetime effect and is chosen to be two Gaus- 
sians convoluted with two double-sided exponentials. Fi- 
nally, the PDF for the remaining events, 'P"('^g^, is the 
convolution of an exponential function with an effective 
lifetime and two Gaussians. 



The fractions (f^hcJif 



isbc^Jidir and /Xr) of sbc, Idlr and 
other events, are determined directly with the B^B'^ and 
B^B~ Monte Carlo simulation. The fractions /"j,'^ of 
obc events are fitted to the data, constraining the v&- 
f obc/ f obc to the estimate obtained with Monte Carlo 

samples. The fraction of B^B~ events is determined 

from the data themselves. 

The last component of the dilepton sample originates 
from non-BB events, mainly continuum events, and has 
been estimated using off-resonance events to represent a 
fraction fcont = (3.1±0.1)% of the data set. To model its 
PDF we use off-resonance events with looser cuts and on- 
resonance events that fail the continuum-rejection cut on 
the Fox- Wolfram moment ratio. The charge asymmetries 
a™^' obtained with the two samples are consistent with 
zero at the 1% level and thus are fixed to zero in the 
likelihood. 



The likelihood is given by 

£(At) = (1 + <Jiaf°"*)(l + q2ai°"^) fcontVcont 

+ (1 - .fcont){f+-VB+B~ + (1 - /+-)^B0S«} 

'PgOgO = fsig){l + qiai'^'"'){l + q2ai^'"')'Pcasc 

I /I I dir\ / -t . dir\^n 

B+B- = (1 - + aiafr°'')(l + a2at2''*'')7'ca3c 

+ 9iafi"')(l + q2afl'')'Plig 

/n,c ^nn,c I r 
other other ' '. 



V, 



case 



IdlT ' IdlT + J sbc ' sbc ~^ J obc ' obc ' 



where qi, q2, fi and f2 are the charges and the flavors 
(e, /i) of the two leptons. 

The likelihood fit gives \q/p\ - 1 = (-0.8 ±2.7) x lO^^, 
Imz = (-13.9 ± 7.3) x lO'^, and AT x Rez = (-7.1 ± 
3.9) X 10^'^ ps^^. The correlation between the measure- 
ments of Imz and AT x Rez is 76 %. If we fix AF = 0, 
we obtain Imz — (—3.7 ±4.6) x 10"'^. Figure^shows the 
At /CP asymmetry between (i?+,£+) and {£^,£~) dilep- 
tons defined in Eq. |21 and the A^prp asymmetry be- 
tween {i^ , £~ ) dileptons with At > and At < defined 
in Eq. 01 




12 14 

^t| (ps) 



FIG. 1; (a) At/cp asymmetry between {£^,£^) and {£~ ,£~). 
A larger charge asymmetry for cascade muons, dominant 
at small |At|, explains the non-flatness of the curve, (b) 
AcpT/CP asymmetry between {£'^ ,£~) dileptons with At > 
and At < 0. 



The T/CP and CPT/CP violation parameters are ex- 
tracted from a binned maximum likelihood fit of the 
events that pass the dilepton selection. The likelihood 
C combines the charge asymmetries in detection and 
the time-dependent PDFs described previously. As the 
charge asymmetries are significantly different for elec- 
trons and muons, we split the sample into four lepton 
combinations: ee, e/x, /Lte and /z/i, in which the first lep- 
ton has the higher momentum. 



There are several sources of systematic uncertainty in 
these measurements. To determine their effect, we vary 
each source of systematic uncertainty by its known or 
estimated uncertainty, and take the resulting deviation 
in the CP parameter as its systematic uncertainty. 

For \q/p\, the most important systematic errors are 
due to uncertainties on electron charge asymmetries. A 
1.4 X 10"'^ deviation of \q/p\ is observed by shifting simul- 
taneously the electron charge asymmetries by 1.0 x 10"'^ 
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which corresponds to the uncertainty estimated with 
Monte Carlo and control samples. The systematic un- 
certainty related to the charge asymmetry due to the 
tracking is estimated by randomly removing a fraction 
equal to 1.6 x 10"'^ of the negative tracks from our data 
sample. This fraction has been determined from an inde- 
pendent data control sample. A 1.0 x 10~^ deviation of 
\q/p\ is observed. Similarly, the 1% uncertainty on charge 
asymmetry for non-BB backgrounds induces a system- 
atic error of 0.6 x 10~^. 

The widths of the first and second Gaussian of the 
resolution function for the obc and sbc categories as well 
as the pseudo-lifetime for the Idlr and other categories 
are varied separately by 10%. This variation is motivated 
by the comparison of the fitted parameters of the signal 
resolution function obtained on generic BB Monte Carlo 
samples and on data being in agreement at 10% level. 
The fractions of the short-lived and long-lived charmed 
meson components for obc and sbc are varied by 10 %. 

We have also varied the parameters Am, t^o and Tg± 
independently within their known uncertainties |lfij | and 
Ar from 10~^ to 0.1. Finally, one of the dominant sys- 
tematic errors on AF x Re z is imperfect knowledge of 
the absolute z scale of the detector and the residual un- 
certainties in the SVT local alignment, giving an error of 
1.2 X 10-3 ps-i. 

TABLE I: Summary of systematic errors for \q/p\, Imz, and 
Ar X Rez measurements. 



Systematic Effects a{\q/p\) cr(Imz) CT(ArxRez) 
(xlO~^) (xlO~^) (xlO'^ps'^) 



Ch. asym. of non-BB bkg 0.6 0.0 0.0 

Ch. asym. in tracking 1.0 0.0 0.0 

Ch. asym. of electrons 1.4 0.0 0.0 

PDF modeling 0.3 2.5 1.2 

Fraction of bkg components 0.2 0.4 0.1 

Am, Tso, Tgi and AF 0.2 1.9 1.1 

SVT alignment 0.5 0.6 L2 

Total Tg 3^2 To 



For each parameter, the total systematic error is the 
sum in quadrature of the estimated systematic errors 
from each source, as summarized in Tabled When we as- 
sume AF = 0, the systematic error for Imz is2.9xl0-'^. 

If we compare our results to AF x Rez = 0.0 and 
Imz = 0.0 (no CPT violation case), the is 3.25 for 
2 degrees of freedom, which gives a confidence level of 
19.7%. Finally, assuming AF = 0, we obtain Imz = 
(-3.7 ± 4.6(stat.) ± 2.9(syst.)) x lO'^. 

In summary with the 1999-2004 data (232 x 10^ BB 
pairs), we have performed a simultaneous likelihood fit of 
the same-sign and opposite-sign dileptons. We measure 
the independent parameters governing CP and T viola- 



tion, and the CPT and CP violation parameters. The 
results are 

\q/p\-l = (-0.8±2.7(stat.)±1.9(syst.)) X 10"^ 

Imz = (-13.9±7.3(stat.)±3.2(syst.)) X 10-^ 

AFxRcz = (-7.1 ±3.9(stat.) ±2.0(syst.)) X lO^^ps-^ 

These measurements are a clear improvement over the 
most precise results previously published Sim. The 
new measurement of \q/p\ is consistent with the Standard 
Model predictions 0. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the computing 
organizations that support BABAR. The collaborating 
institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DEC (Germany), 
INFN (Italy), EDM (The Netherlands), NER (Norway), 
MIST (Russia), and PPARC (United Kingdom). Indi- 
viduals have received support from CONACyT (Mex- 
ico), Marie Curie EIF (European Union), the A. P. Sloan 
Foundation, the Research Corporation, and the Alexan- 
der von Humboldt Foundation. 



* Also at Laboratoire de Physique Corpusculaire, 
Clermont-Ferrand, France 

t Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

* Also with Universita della Basilicata, Potenza, Italy 

[1] J. H. Christenson, J. W. Cronin, V. L. Fitch and R. 

Turlay, Phys. Rev. Lett. 13, 138 (1964). 
[2] CPLEAR Collaboration, A. Apostolakis et al, Phys. 

Lett. B 456, 297 (1999). 
[3] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 70, 

012007 (2004). 

[4] M. Beneke et al, Phys. Lett. B 576, 173 (2003); M. 

Ciuchini et al, JHEP 0308, 031 (2003). 
[5] BABAR Collaboration, B. Aubert et al, Nucl. Instr. 

Methods Phys. Res., Sect. A 479, 1 (2002). 
[6] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

88, 231801 (2002). 
[7] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 

(1978). 

[8] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

88, 221803 (2002). 
[9] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 66, 

032003 (2002). 

[10] Heavy Flavor Averaging Group, K. Anikeev et al, 

hep-ex/0505100 
[11] BELLE Collaboration, E. Nakano et al, submitted to 

Phys. Rev. D , hep-ex/0505017. 



